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Andreev current enhancement and subgap conductance of superconducting hybrid structures in the 

presence of a small spin-splitting field 
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We investigate the subgap transport properties of a S-F-Ne structure. Here S (Ne) is a superconducting 
(normal) electrode, and F is either a ferromagnet or a normal wire in the presence of an exchange or a spin- 
splitting Zeeman field respectively. By solving the quasiclassical equations we first analyze the behavior of 
the subgap current, known as the Andreev current, as a function of the field strength for different values of the 
voltage, temperature and length of the junction. We show that there is a critical value of the bias voltage V * above 
which the Andreev current is enhanced by the spin- splitting field. This unexpected behavior can be explained 
as the competition between two-particle tunneling processes and decoherence mechanisms originated from the 
temperature, voltage and exchange field respectively. We also show that at finite temperature the Andreev 
current has a peak for values of the exchange field close to the superconducting gap. Finally, we compute the 
differential conductance and show that its measurement can be used as an accurate way of determining the 
strength of spin- splitting fields smaller than the superconducting gap. 

PACS numbers: 74.45.+c,74.25.F- 



Introduction- Transport properties of hybrid structures con- 
sisting of superconducting and non- superconducting materials 
have been studied extensively in the last decades ^ In partic- 
ular, there is a renewed interest in the study of the subgap con- 
ductance of superconducting hybrid structures in the presence 
of Zeeman-like fields in view of the presumable detection of 
Majorana Fermions^. Intuitively, due to the gap A in the den- 
sity of states of a superconductor the charge transport through 
a superconductor-normal (S-N) metal junction is expected to 
vanish for voltages smaller than A. However, this is not al- 
ways the case. Experiments on S/N structures have shown 
a finite subgap conductance^. This behavior was discussed 
theoretically in Refs. 4 and 5. It was shown that the conduc- 
tance of a S-N-Ne structure, where Ng denotes a normal metal 
electrode, shows a peak at a voltage smaller than the super- 
conducting gap^ in the case of finite S/N barrier resistances or 
if N is a diffusive metal^'^ . A similar behavior was predicted 
if one substitutes the normal by a ferromagnetic metal (F)2r-. 
In all these examples, the key mechanism to explain the finite 
subgap conductance is the Andreev reflection ^^-^^ which takes 
place at the S/N and S/F interfaces and allows the flow of an 
electric current even for voltages smaller than the supercon- 
ducting gap A. By this process an electron from the normal 
region is reflected as a hole forming a coherent electron-hole 
pair which penetrates into a diffusive normal region over dis- 
tances of the order of the thermal length ^/Wff , where Q) is 
the diffusion coefficient and T is the temperature (here and 
below we set /i = ^5 = 1). This mechanism leads to a finite 
condensate density in the normal metal, i.e. to the so called 
superconducting proximity effect. 

At a S/F interface the mechanism of charge transport is 
however modified since the incoming electron and reflected 
hole belong to different spin bands ^^. Thus, one expects a sup- 
pression of the Andreev current by increasing the exchange 



field h of the ferromagnet, which is a measure for the spin- 
splitting at the Fermi level. In the ferromagnet the coherence 
length of the electron-hole pairs is given by the m inimum 
between the thermal and the magnetic (~ ^/WJh) lengths. 
One expects that by increasing the strength of the field h the 
electron-hole coherence would be suppressed and hence the 
subgap current reduced. As we show below, this intuitive pic- 
ture does not hold always. 

In this Letter we analyze the Andreev current and conduc- 
tance through a S-F-Ne hybrid structure as a function of the 
field h. Here h denotes either the intrinsic exchange field of a 
ferromagnet or a spin- splitting field in a normal metal caused 
by either an external magnetic field or the proximity of an in- 
sulating ferromagnet^^. We focus our study on weak fields, 
/z ^ A and h> A. We find an interesting interplay between 
phase-coherent diffusive propagation of Andreev pairs due to 
the proximity effect and decoherence mechanisms originated 
from the temperature, voltage and exchange field respectively. 
This interplay leads to a non-monotonic behavior of the trans- 
port properties as a function of h. For very low temperatures 
and voltages eV <^ A the Andreev current decays monotoni- 
cally by increasing h as expected. If one keeps the voltage low 
but now increases the temperature, the Andreev current shows 
a peak at /z ^ A. An unexpected behavior is obtained when 
the voltage exceeds some critical value V*. In this case, the 
Andreev current is enhanced by the field h reaching a max- 
imum 8it h ^ eV. We show that the value of V* depends on 
the length of the F wire and the temperature. In particular, for 
zero-temperature and in the long-junction limit, i.e. when the 
length of F is much larger than the coherence length, we show 
that eV* ^ 0.56Ao, where Aq is the value of A at T = 0. We 
also compute the subgap conductance of the system at low 
temperatures and small fields h < A. We show that it has a 
peak at eV = h. Thus, transport measurements of this type 



can be used to determine the strength of a weak exchange or 
Zeeman-Hke field in the nanostructure. 

Model and basic equations- We consider a ferromagnetic 
wire F. Its length, L is smaller than the inelastic relaxation 
length. The wire is attached at x = to a superconducting (S) 
and at X = L to a normal (Ne) electrode. As noticed above, F 
can also describe a normal wire in a spin- splitting field B (in 
which case h = {IbB, where /i^ is the Bohr magneton) or in 
proximity with an insulating ferromagnet^^. We consider the 
diffusive limit, i.e. we assume that the elastic scattering length 
is much smaller than the decay length of the superconducting 
condensate into the F region. In order to describe the transport 
properties of the system we compute the quasiclassical Green 
f unctions ^^'^^. They obey the Usadel equation^^ that in the so 
called 6 — parametrization reads^^ 



5xx^± = 2/ 



E±h 



- sinh 6± . 



(1) 



Here the upper (lower) index denotes the spin-up (down) com- 
ponent. The normal and anomalous Green functions are given 
by g± = cosh 6± and f± = sinh 6± respectively. Because of 
the high transparency of the F/Ne interface the functions G± 
vanish at x = L, i.e. superconducting correlations are negligi- 
ble at the F/Ne interface. We consider a tunneling barrier at 
the S/F interface and assume that its tunneling resistance Rj 
is much larger than the normal resistance Rf of the F layer. 
Thus, by voltage-biasing the Ng the voltage drop takes place 
at the S/F tunnel interface. To leading order in Rp/Rj <C 1 
the Green functions obey the Kupriyanov-Lukichev boundary 
condition at x = 0^^ 



dxG±\x= 
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where 6s = arctanh(A/£^) is the superconducting bulk value 
of the function 6. Once the functions 6± are obtained one 
can compute the current through the junction. In particular, 
we are interested in the Andreev current, i.e. the current for 
voltages smaller than the superconducting gap due to Andreev 
processes at the S/F interface. Such current is given by the 
expression 
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(3) 
where n-{E) = ^(tanh[(£ + eV)/2T] - tanh[(£ - eV)/2T]) 
is the quasiparticle distribution function in the Ng electrode, 
a±{E) = {l/^)J^dx cosh-^[Ree±(x)], W = ^Rf/2LRt is 
the diffusive tunneling parameter ^^'^^, and ^ = y/Wj2A is the 
superconducting coherence length. Eq. ^ is the expression 
used throughout this article in order to determine the subgap 
charge transport—. 

Results- We first compute the Andreev current numerically 
by solving Eqs. ([T]|3]). In Fig. [T] we show the dependence of the 
Andreev current on the exchange field h for different values of 
the bias voltage and temperature for a ferromagnetic F wire of 
the length L= 10(^. 

We consider first the zero-temperature limit. For small 
enough voltages {e.g eV = 0.3 A, black solid line in Fig. (Tt) 
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FIG. 1. (Color online) The /^-dependence of the ratio /^(/?)//^(0) for 
L=\O^SindW = 0.007. Left panels correspond to (a) eV = 0.8A and 
(b) eV = 0.3 A. The different curves are for T = (black solid line), 
T = 0.12Ao (blue dashed line) and T = 0.25Ao (red point-dashed 
line). The right panels corresponds to (c) T = 0.25Ao and (d) T = 0, 
while the different curves to eV = 0.3 A (black solid line), eV = 0.55 A 
(blue dashed line) and eV = 0.8A (red point-dashed line). In the 
(a) panel curves are vertically shifted with respect to each other for 
clarity. 



the Andreev current decays monotonously with increasing h. 
This behavior is the one expected, since by increasing h the 
coherence length of the Andreev pairs in the normal region is 
suppressed, leading to a reduction of the subgap current. For 
large enough voltages (e.g eV = 0.8A in Fig.[T]) and keeping 
the temperature low, the Andreev current first increases by in- 
creasing h, reaches a maximum at h ^ eV, and then drops by 
further increase of the exchange field, as it is shown for exam- 
ple by the black solid line in Fig.fT^. A common feature of all 
the low- temperature curves in Fig. [T]is the sharp suppression 
of the Andreev current ath^eV. 

For large enough temperatures (T = 0.25Ao in Fig[T]:) one 
observes a peak at h ^ A [Fig. [He)]. The relative height of 
this peak increases with temperature and voltage as one sees 
in Figs. [It and[T]: respectively. In the case of large enough 
values of V and T, one is able to observe both the enhance- 
ment of the Andreev current by increasing h and the peak at 
h^ A (see for example blue dashed line in Fig[Tt). For val- 
ues of the exchange field larger than A, the Andreev current 
decreases by increasing h in all cases . In principle all the 
behaviors of the Andreev current can be observed by measur- 
ing the full electric current through the junction as the single 
particle current is almost independent of h. 

In order to give a physical interpretation of these results, we 
first recall the details of the process of two-electron tunneling 
that gives rise to subgap current"^ in diffusive systems in the 
absence of an exchange field. The value of this current is gov- 
erned by two competing effects. On the one hand, the origin 
of the subgap current is the tunneling from the normal metal to 
the superconductor of two electrons with energies ^^^ and ^^^ , 
respectively and momenta ki and k2, that form a Cooper pair. 
This process is of the second order in tunneling and there- 
fore involves a virtual state with an excitation on both sides 
of the tunnel barrier. The relevant virtual state energies are 

given by the difference E^ — ^kx ,k2 ' where Ej^ = JA^ + ^^ is 
the excitation energy of a quasiparticle with the momentum 
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FIG. 2. (Color online) Schematic energy diagram for a non magnetic 
(a) and magnetic (b) metal. The thick solid parabolas are the disper- 
sions of free electrons with spin up (t) and spin-down (|). The k axis 
corresponds to the Fermi level in the superconductor. We consider 
quasi-electrons and -holes with energies i^V < A and momentum 
^12 . Time-reversal pairing requires that ki =k2. In case of a normal 
metal [panel (a)] this condition is satisfied only for eV = while for 
a ferromagnet (hy^O)ifh = eV. 



k in the superconductor. Typical values of ^^ ^^^ T or eV. 
Hence under subgap conditions T, eV <C A, the virtual state 
energy is typically given by the superconducting gap A. How- 
ever, when these characteristic energies become larger and ap- 
proach the value of the gap, the difference E^ — ^kiM even- 
tually vanishes. As a result, the amplitude for two-electron 
tunneling increases drastically, leading to a strong increase 
of the Andreev current, accompanied by the onset of single- 
particle tunneling at energies above the gap A. On the other 
hand, two-electron tunneling is a coherent process: the main 
contribution to two-electron tunneling stems from two nearly 
time-reversed electrons ki c^ —k2 located in an energy win- 
dow of width 5£ ~ eV^ T close to the Fermi energy, diffusing 
phase-coherently over a typical distance Lcoh = \/^/(5£) in 
the normal metal before tunneling^. This coherence length de- 
creases upon increasing the characteristic energies Sec^eV^T, 
thereby decreasing the Andreev current. 

We now turn to the effect of the exchange field h on two- 
electron tunneling. If h is nonzero, the majority and minority 
spin electrons at the Fermi level are characterized by differ- 
ent wave vectors kf^± = kf ^ 5k, where 5^ ~ h/vf and vf 
is the Fermi velocity. In Fig. [2] we show a schematic energy 
diagram. The wave vectors ^77 ± are determined by intersec- 
tion between the parabolas and the /:-axis. For a given value 
of ^y < A and in the absence of an exchange field the rele- 
vant excitations with energies ~ ±eV and wave vectors ^1^2 
are not time-reversed (see Fig. [2^) and therefore their contri- 
bution to the current is not coherent. However upon increas- 
ing h, \ki\ -^ \k2\, i.e the relevant excitations become more 
and more coherent, leading to an additional increase of two- 
electron tunneling. In particular when h = eV , ki = —k2 (cf. 
Fig. [2I3). If r ^ there are no occupied states for ^]^> eV. 
Consequently as soon sls h> eV, the energy window around 
the Fermi level does not contain time-reversed electrons. This 
leads to the drop of the Andreev current shown for example in 
Fig. [T]l. In contrast, for finite values of T, there are thermally 
induced quasiparticles with energy ~ A, that become exactly 
time-reversed whenever h = A. This leads to the maximum of 
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FIG. 3. (Color online) Voltage-junction length (a) and voltage- 
temperature (b) diagrams. The black solid line represents the values 
of eV* /A. For the range of parameters situated below this line the 
Andreev current decreases in the presence of a small exchange field 
(suppression), while in the region above the line the current increases 
(enhancement). We set W = 0.007 in both panels, T = in panel (a) 
and L= \0^ in panel (b). 

the current at h = A when the temperature is finite (cf. Fig. 
[T]:). The effects are most clearly seen when plotting the ra- 
tio lA{h)/lA{0), as the Andreev pair decoherence effects due 
to temperature or voltage are then divided out. 

A more quantitative understanding of the effects discussed 
above can be get by analyzing some limiting cases in which 
simple analytical expressions for the current can be derived. 
We first focus our analysis on the zero-temperature limit. Due 
to the tunneling barrier at the S/F interface the proximity effect 
is weak and hence one can linearize Eqs. ([TIO with respect to 
Rf/Rt <C 1. After a straightforward calculation one obtains 
the Andreev current in this limit, 
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For a large exchange field, h:^ Aq> eV one can evaluate this 
expression obtaining 
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(5) 



Thus, the Andreev current decays as h~^^^ for large values of 
h in accordance with our numerical results (see Fig. [T). 

In the case of small values of /z, /z < ^V < Aq, one can eval- 
uate Eq. dU) in the long-junction limit, i.e. when L ^ \/D/h. 
In this case the Andreev current reads 
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Aq^Rf 
eLR\ 



arctanh f 



eV±jh\ 
Ao+ jh J 
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This expression describes the two different behaviors obtained 
in Fig. \l\for h < eV. For small voltages I a decreases by 
increasing the field h. However, for large enough values 
of the voltage I a is enhanced by the presence of the field. 
From Eq. (O we can determine the voltage V*, at which the 




FIG. 4. (Color online) The bias voltage dependence of differential 
conductance at T = for fields: h = 0.3 (black solid) and h = 0.^ 
(blue dashed). Here Ga = 4RtGa, W = 0.007 and L = lOc^. 

crossover between these two behaviors takes place, by ex- 
panding the expression for the current up to second order in 
h/eV <C 1, i.e. up to the first non- vanishing correction to Ia 
due to the exchange field. This expansion leads to the follow- 
ing transcendental expression which determine the voltage V * 
at which the crossover takes place, 

^ ) ^2 (^^^^^^^ \/^^*Mo + arctan V^y*/Ao) . 

(7) 
From here we get eV ^ 0.56Ao. For V < V* the Andreev cur- 
rent decays monotonically with h while for V > V * it increases 
up to a maximum value 3.th<eV. This is in agreement with 
our numerical results in Fig. [T] 

For an arbitrary length L and finite temperature we have 
computed the value of V* numerically. In Fig. [3] we show the 
results. The solid black line gives the values of V* as a func- 
tion of L and T [the (a) and (b) panels of Fig. [3] respectively]. 
The area below the black curve corresponds to the range of 
parameters for which the Andreev current is suppresses by 
the presence of a spin- splitting field, while the area above the 
solid line corresponds to the range of parameters for which the 
unexpected enhancement of the subgap current takes place. 
According to Fig.[3a) at T = the value of V* first decreases 
as L increases, reach a minimum and then grows again up to 
the asymptotic value ^V* ^ 0.56Ao. Also the dependence of 
y* on the temperature is non-monotonic having a maximum 
value at r ~ 0.2Ao. 

Small spin- splitting fields, as those studied in the present 
work, can be created by applying an external magnetic field B, 
in which case h = jIbB or by the proximity of a ferromagnetic 
insulator as discussed in Ref . [Tsl It may be also an intrinsic 
exchange field of weak ferromagnetic alloys (see, for exam- 
ple, Ref. 2JJ. Such small exchange fields are in principle dif- 



ficult to detect. However, as we show in Fig. HI by measuring 
the subgap differential conductance G = dl/dV at low temper- 
atures, one can accurately determine the value of h. AtT = 
the conductance shows two well defined peaks, one aieV = h 
and the other at eV = A. These are related to a sudden increase 
of the coherence length between the electron-hole pairs in the 
ferromagnet and of the two-particle tunneling amplitude re- 
spectively. As we have seen above, at small voltages eV < h 
electrons with majority spins do not find time-reversed part- 
ners in the narrow energy window around the Fermi energy, 
i.e. such pairs show weak coherence. By increasing the volt- 
age eV <h, the contribution of time-reversed electrons to the 
current gradually increases and consequently the differential 
conductance increases, reaching a maximum at eV = h. Fur- 
ther increase of the voltage, eV > h, leads to an increasing 
contribution to the current from non time-reversed electron- 
hole pairs and therefore to a suppression of the coherent con- 
tribution to G. At /z < ^y < A the two-electron tunneling am- 
plitude increases as {eV — A)~^ due to virtual state contribu- 
tions with energies eV close to the gap; as a result the con- 
ductance shows a sharp increase. For h ^ (normal metal) 
the peak moves toward eV ^0 (not shown here), which cor- 
responds to the zero bias peak discussed, for example, in Ref. 

In conclusion, we present an exhaustive study of the sub- 
gap charge current through S-F-Ne hybrid structures in the 
presence of a spin- splitting field. We have demonstrated the 
existence of a threshold bias voltage y* above which the An- 
dreev current can be enhanced by an exchange field. We also 
have shown that at finite temperatures the Andreev current has 
a peak for values of the exchange field close to the supercon- 
ducting gap A. Finally, we have proposed a way to determine 
the strength of small exchange fields by measuring the dif- 
ferential conductance. Beyond the fundamental interest, our 
results can also be useful for the implementation of a recent 
and interesting proposal ^^ which suggests a way to detect the 
odd-triplet component^^ of the superconducting condensate 
induced in a normal metal in contact with a superconductor 
and a ferromagnetic insulator. The latter induces an effective 
exchange field in the normal region. The amplitude of such 
induced exchange fields is smaller than the superconducting 
gap2i. Therefore the proposed ferromagnet proximity system 
in Ref. 43' is a candidate to observe the phenomena described 
in the present work. 
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